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Abstract 

New physical and mathematical applications of recently invented 
fractional Poisson probability distribution have been presented. 

As a physical application, a new family of quantum coherent states 
have been introduced and studied. 

Mathematical applications are related to the number theory. We 
have developed fractional generalization of the Bell polynomials, the 
Bell numbers, and the Stirling numbers of the second kind. 

The fractional Bell polynomials appearance is natural if one evalu- 
ates the diagonal matrix element of the evolution operator in the basis 
of newly introduced quantum coherent states. 

The fractional Stirling numbers of the second kind have been intro- 
duced and applied to evaluate skewness and kurtosis of the fractional 
Poisson probability distribution function. A new representation of 
the Bernoulli numbers in terms of fractional Stirling numbers of the 
second kind has been found. 

In the limit case when the fractional Poisson distribution becomes 
the well-known Poisson probability distribution all of the above listed 
new developments and implementations turn into the well-known re- 
sults of the quantum optics and the number theory. 



* E-mail address: nlaskin@rocketmail.com 



1 



PACS numbers: 05.10.Gg; 05.45.Df; 42.50.-p. 

Keywords: fractional Poisson process, quantum coherent states, 
generating function, generalized Stirling and Bell numbers. 

1 Introduction 

In the past decade it has been realized that the understanding of complex 
quantum and classical physic phenomena has required the implementation 
of long-range space and long memory processes. 

The mathematical model to capture the long-range space dependence im- 
pact on quantum phenomena is the Levy path integral approach invented and 
studied in Ref. PQ, [2]. The Levy path integrals generalize the path integral 
formulation of the quantum mechanics developed in 1948 by Feynman [3], 

0- 

To study a long memory impact on the counting process, the fractional 
Poisson random process had been introduced and developed in [5]. The 
fractional Poisson distribution captures the long-memory effect which re- 
sults in non-exponential waiting time probability distribution function em- 
pirically observed in complex quantum and classical systems. The quantum 
system example is the fluorescence intermittency of single CdSe quantum 
dots, that is, the fluorescence emission of single nanocrystals exhibits in- 
termittent behavior, namely, a sequence of "light on" and "light off' states 
departing from Poisson statistics. In fact, the waiting time distribution in 
both states is non-exponential [B]. As examples of classical systems let's men- 
tion the distribution of waiting times between two consecutive transactions 
in financial markets [7] and another, which comes from network communica- 
tion systems, where the duration of network sessions or connections exhibits 
non-exponential behavior [8]. 

The non-exponential waiting time distribution function had been ob- 
tained for the first time in [PJ, based on the fractional generalization of the 
Poisson exponential waiting time distribution. 

The fractional Poisson process is a generalization of the well known Pois- 
son process. A simple analytical formula for the fractional Poisson probability 
distribution function has been obtained in [S] based on the fractional gen- 
eralization of the Kolmogorov- Feller equation introduced in [5J. Comparing 
to the standard Poisson distribution, the probability distribution function of 
the fractional Poisson process [5] has an additional parameter /i, < /i < 1. 
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In the limit case fi = 1 the fractional Poisson becomes the standard Pois- 
son process and all our findings are transformed into the well-known results 
related to the standard Poisson probability distribution. 

Now we present quantum physics and the number theory applications of 
the fractional Poisson probability distribution. 

The quantum physics application is an introduction of a new family of 
quantum coherent states. The motivation to introduce and explore the new 
coherent states is the observation that the squared modulus | < n\q > | 2 of 
projection of the newly invented coherent state |s > onto the eigenstate of 
the photon number operator |n >, gives us the fractional Poisson probability 
Pfj,{n) that n photons will be found in the new coherent state \<; > . Following 
the Klauder's framework to qualify quantum states as generalized coherent 
states [10] , we prove that our quantum coherent states, (i) are parametrized 
continuously and normalized; (ii) admit a resolution of unity with positive 
weight function; (iii) provide temporal stability, that is, the time evolution 
of coherent states remains within the family of coherent states. We have 
defined the inner product of two vectors in terms of their coherent state |^ > 
representation and introduced the functional Hilbert space. 

Mathematical applications are related to the number theory. The Bell 
polynomials, the Bell numbers [11] and the Stirling numbers of the second 
kind [12] - [TJ] have been generalized based on the fractional Poisson proba- 
bility distribution. 

Appearance of the fractional Bell polynomials is natural if one evaluates 
the diagonal matrix element of the quantum evolution operator in the basis 
of newly introduced quantum coherent states. The appearance of fractional 
Bell numbers manifests itself in the fractional generalization of the celebrated 
Dobinski formula [T5|, [T5] for the generating function of the Bell numbers. 

Fractional Stirling numbers of the second kind have been applied to eval- 
uate skewness and kurtosis of the fractional Poisson probability distribution. 

In the limit case when \i = 1, the fractional Poisson distribution becomes 
the well-known Poisson probability distribution, and all above listed new 
developments and findings turn into the well-known results of the quantum 
coherent states theory [T7], [TS], [IS] and the number theory [IB] , [Tij . 

The paper is organized as follows. 

The basic definitions of the fractional Poisson random process are briefly 
reviewed in Sec. 2, where Table 1 has been presented to compare formulas 
related to the fractional Poisson process [5] to those of the well-known ones 
related to the standard Poisson process. In Sec. 3 we introduce and study new 
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quantum coherent states and their applications. Fractional generalizations 
of the Bell polynomials, Bell numbers and Stirling numbers of the second 
kind have been introduced and developed in Sec. 4. New equations for the 
generating functions of fractional Bell polynomials, fractional Bell numbers 
and fractional Stirling numbers of the second kind have been introduced and 
elaborated. The relationship between the Bernoulli numbers and fractional 
Stirling numbers of the second kind has been found. Table 3 presents a few 
fractional Stirling numbers of the second kind. Statistics of the fractional 
Poisson probability distribution were studied in Sec. 5. The centered moment 
of m-order has been obtained in terms of fractional Stirling numbers of the 
second kind. Variance, skewness and kurtosis of the fractional Poisson prob- 
ability distribution function have been presented in terms of the centered 
moments of m-order. We also explain and discuss how the well-know equa- 
tions of the quantum optics and the number theory, related to the standard 
Poisson probability distribution, can be obtained from our generalized re- 
sults. Tables 2 and 4 summarize the new results related to the fractional 
Poisson distribution vs those related to the standard Poisson probability dis- 
tribution. 

2 Fundamentals of the fractional Poisson pro- 
cess 

The fractional Poisson process has been introduced and developed [5] as the 
counting process with probability P M (n, t) of arriving n items (n = 0, 1, 2, ...) 
by time t. The probability P^(n,t) is governed by the fractional generaliza- 
tion of the Kolmogorov-Feller equation 

L>?P>, t) = u (P> - 1, t) - P„(n, t)) + 6 n>0 , < fi < 1, (1) 

with normalization condition 

oo 

X>„(M) = 1, (2) 

n=0 

where qD^ is the operator of time derivative of fractional order fx defined 
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as the Riemann-Liouville fractional integral, 

drf(r) 



T(-fi) J (t-r)^' 
o 

and S nt Q is the Kronecker symbol, the gamma function T(n) has the fa- 

oo 

miliar representation r(/x) = j dte^ t t fJ- ^ 1 , Re/x > 0, and parameter v has 

o 

physical dimension [u] = sec~ M . The initial condition P^(n,t = 0) = <5 n .o is 
incorporated in Eq.(pQ). One can consider the fractional differential-difference 
equation ([1]) as a generalization of the differential-difference equation which 
defines the standard Poisson process (see, for instance, Eq (1.4.4) in Ref. |23j). 

To solve Eq.([T]) it is convenient to use the method of the generating 
function. We introduce the generating function G M (s,t) 

oo 

G>,t) = £VP M (n,*). (3) 

n=0 

Then by multiplying Eq.([T]) by s n , summing over n, we obtain the follow- 
ing fractional differential equation for the generating function G^s, t) 



v n=0 n=0 

z/(s-l)G M (s,t)+ / ' 



r(i-M) 

The solution of this fractional equation has a form 

G li (s,t)=E ll (vtP(s-l)), (5) 

where E^z) is the Mittag-Leffler function given by its power series [24J, 
[250 

m=0 V ^ y 



1 The basic formulas on fractional calculus can be found in Refs. [20] - (22] . 
2 At /x = 1 the function E^(z) turns into exp(z). 
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Expanding (jSJ) in series over s results in accordance with definition ([3]) 



/ s v 2 ^ (k + n)\ {-vt^) k 

P ^*)-hr^ k\ T{p(k + n) + 1)) ' 0< ^ L (7) 

The P M (n, i) gives us the probability that in the time interval [0, t] we 
observe n counting events. When \i = 1 the P M (n, t) is transformed to the 
standard Poisson distribution (see Eq.(14) in Ref.[5]). Thus, Eq.([7j) can be 
considered as a fractional generalization of the standard Poisson probability 
distribution. The presence of an additional parameter // brings new features 
in comparison with the standard Poisson distribution. 

On a final note, the probability distribution of the fractional Poisson 
process can be represented in terms of the Mittag-Leffler function E^(z) in 
the following compact way [5], 

(—z) n d n 

P,(n,t) = y —^—E,(z)\ z= ^ (8) 

P lt (n = 0,t) = £„(-!/#*). (9) 

At /j, = 1 Eqs.flEJ) and ([9]) are transformed into the well known equations 
for the standard Poisson process with the substitution v — > 77, where 77 is 
the rate of arrivals of the standard Poisson process with physical dimension 
77 = sec -1 . 

Table 1 summarizes equations attributed to the fractional Poisson process 
with those belonging to the well-known standard Poisson process. Table 1 
presents two sets of equations for the probability distribution function P(n, t) 
of n events having arrived by time t, the probability P(0, t) of nothing having 
arrived by time t, mean n, variance a 2 , generating function G(s,t) for the 
probability distribution function, moment generating function H(s,t), and 
waiting-time probability distribution function ip{r). 
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fractional Poisson (0 < /i < 1) 


roisson [fi = 1) 


P(n,t) 


oo 

(vt*) n (k+n)\ 




( J exp( vt) 


nl ^ k\ r(/Lt(fc+n)+l)) 
fc=0 




P(n,t) 




\z=—vt» 


{ J ^»exp(«)lz=-^ 


P(0,t) 




exp(— z/t) 


n 




z/t 


a 2 


i/f , ( \ 2 

roi+i) 1 ^r^+i),/ 1 






G(s,t) 




-1)) 


exp{z7t(s — 1)} 


H(s,t) 


E^{e- S 


-1)) 


exp{z7t(e _s — 1)} 


V(r) 






Ve~ VT 



Table 1. Fractional Poisson process vs the Poisson proces, 



3 New family of coherent states 

The quantum mechanical states first introduced by Schrodinger [26] to study 
the quantum harmonic oscillator are now well-known as the coherent states. 
Coherent states provide an important theoretical paradigm to study electro- 
magnetic field coherence and quantum optics phenomena [17], |18j . 

The standard coherent states are given for all complex numbers zeC, 

by 



\z >= e 



( 2 a+- 2 *a)| >= e -I N 2 ^ ^—\U >, (10) 



n=0 v "' 



where a + and a are photon field creation and annihilation operators that 
satisfy the Bose-Einstein commutation relation [a, a + ] = aa + — a + a = 1, and 
the orthonormal vector \n >— "^=f( a+ )™|0 >> which is an eigenvector of the 
photon number operator iV = a + a, N\n >= n\n >, < n\n' >= 5 n ^ n >. 

The projection of the coherent state \z > onto the state |n > is 



<n\z>=^=e zl 2 ! 2 . (11) 



3 All definitions and equations related to the fractional Poisson process are taken from 

0. 
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Then the squared modulus of < n\z > gives us the probability P{n) that 
n photons will be found in the coherent state \z >. Thus, we come to the 
well- know result for probability P[n) 



\2n 

• |2 



P(n) = I < n\z > I 2 = ^-e" |z| , (12) 

n\ 

which is recognized as a Poisson probability distribution with a mean 
value \z\ 2 . The value \z\ 2 is in fact the mean number of photons when the 
state is a coherent state \z > 

oo 

\z\ 2 = ^nP(n) =< z\a + a\z > . (13) 

n=0 

We introduce a new family of coherent states \q > 

k >= E { ^^{E^\-^)) 1/2 \n >, (14) 

n=0 

and adjoint states < ^| 

< ?i = E < »l^^(^-Mkl*)) 1/a , (is) 

n=0 v 

where 

E^{-^)^—E,{x)\ x= ^ (16) 

and Epix) is the Mittag-Leffler function defined by Eq.(jHJ), complex num- 
ber q stands for labelling the new coherent states, and the orthonormal vector 
| n > is the same as for Eq. ffTUl) . 

To motivate the introduction of new coherent states |<j > we calculate the 
projection of the coherent state \q > onto the state |n >, 

< Ilk >=<4^( £ « ( - Mk i 2 ")) 1/2 , 



then the squared modulus of < n|<j > gives us the probability P^{n) that 
n photons will be found in the quantum coherent state \q >. Thus, we come 
to the fractional Poisson probability distribution of photon numbers P^n) 

P» = | < n| f > I 2 = (E^(-^n) , (17) 
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with mean value (/x|^| 2/i ) /r(/x + 1). The value (/i|^| 2M )/r(/i + 1) is in fact 
the mean number of photons when the quantum state is the coherent state 

k > 

oo 

{fi\^)/r^ + 1) = nP M) =< ^l a+a k > • 

ra=0 

It is easy to see that when /x = 1 we have E\{x) = exp(— x) and Ej?\—x) = 
exp(— x). Hence, Eq. (fllj) turns into to Eg. ([TDD an d Eq. (TTTT) leads to 
Eg. ffTTl) . In other words, the new coherent states defined by Eg. (1141) gener- 
alize the standard coherent states Eg. ffTUl) . and this generalization has been 
implemented with the help of the fractional Poisson probability distribution. 
The new family of coherent states Eg. (|14[) has been designed here to study 
physical phenomena where the distribution of photon numbers is governed 
by the fractional Poisson distribution Eg. (1 171) . 

Let's answer the guestion if the newly introduced coherent states \q > are 
really generalized coherent states? 

Quantum mechanical states are generalized coherent states if they |lUj : 

(i) are parametrized continuously and normalized; 

(ii) admit a resolution of unity with a positive weight function; 

(iii) provide temporal stability, that is, the time evolving coherent state 
belongs to the family of coherent states. 

Let's now show that the new coherent states introduced by Eg. ([T4"|) satisfy 
all above listed conditions. 

To prove (i), we note that the coherent states \q > introduced by Eg. (1141) 
are evidently parametrized continuously by their label q which is a complex 
number q = £ + ir], with £ = Re<^ and r] = Iraq. Because of the normalization 

oo 

condition of the fractional Poisson probability distribution P^{ n ) = 1 an d 

n=0 

< n\n' >= 5 n y, we have 

< qk > = 



EE 

n=0 n'=0 



< n 



*fj,\n 



2^1/2 (v^) ra ^'\-^yf 2 \n' >-. 



;is) 



9 



n=0 ' n=0 

that is, the coherent states \<z > are normalized. 

To prove (ii), that is, the coherent states |^ > admit a resolution of unity 
with a positive weight function, we introduce a positive function W M (|<r| 2 ) > 
which obeys the equation 

1 J d\\,>W,(\,\ 2 )<,\ = I, (19) 

c 

where d 2 q = d(Req)d(lmq) and the integration extends over the entire 
complex plane C. This equation with yet unknown function W-^d^ 2 ) can be 
considered as "a resolution of unity". To find the function W M (|?| 2 ) let's 
transform Eq. (fl9~l) . Introducing of new integration variables p and <fi by 
<; = pe^, d 2 s = p 2 dpd(f> and making use of Eqs. ffT4l) and (jT5l) yield 



i j d\\q > WM 2 ) < q\ 



oo 2tt 

-jrjtf dpp {n+m ^ +1 j d^-^^fl (Ej?\-p\p\ 2 ^) 1/2 x (20) 

n=0 m=0 q q v 

(E^\-p\p\ 2 »)) 1/2 \n><m\ = I. 

By interchanging the orders of summation and integration and carrying 
out the integration over 0, we get a factor 2Ti5 n , m , which reduces the double 
summation to a single one. Therefore, Eq. (l20l is simplified to 

i J d\\q > WM 2 ) < q\ = 

c 

oo 

E "I / dpp 2 ^ +1 W,(p 2 )E^(-p\p\ 2 n\n X n\ = I. (21) 
n=0 U - { 

Because of the completeness of orthonormal vectors |n > 
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CO 



^2\n><n\ = I, (22) 



n=0 



we come to the following integral equation to find the positive function 



00 

J dxx^W^E^i-fix 11 ) = nl (23) 







To solve Eq. (!23|) we use the Laplace transform of the function t^Ejl (—fit ,t ), 
see Appendix, 



f dte-'H^E^i-^) = , n! - S \ - . (24) 
J M (s^ + /i) n+1 v ' 



By comparing Eqs. (l23l) and ( 1241) we conclude that the positive function 



Wfiix) has the form 



W^x) = (1 -Ai)^ -exp{-(l -fi) 1/f *x}, 0</i<l. (25) 

Thus, we proved that the coherent states \q > admit a resolution of unity 
with the positive weight function W^{x) given by Eq. fl25l) . At \i = 1, function 
Wfj,{x) becomes 



Wi(x) = limW / M (x) = lim 



[1 - ji) V • exp{-(l - /i) 1 ^} 



and we come back to the resolution of unity equation for the standard 
coherent states \z > 



d 2 z\z >< z\ = I. (26) 



To prove (iii), we note that if \n > is an eigenvector of the Hamiltonian 
H = hjjN = hu)a + a, where h is Planck's constant, then the time evolution 
operator exp(—iHt/h) results 

exp(-iHt/H)\n >= e~ iujnt \n > . 
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In other words, the time evolution of \n > results in appearance of the 
phase factor only while the state does not change. Let's consider time evo- 
lution of the coherent state \q > defined by Eq. (TH|) . Well, as far as the 
coherent state is not an eigenstate of H then one may expect that it evolves 
into other states in time. However, it follows that 



exp(-iHt/h)\<; >= Y {E^\-^)fl 2 e-^ nt \n >= je"^ >, 



(27) 

which is just another coherent state belonging to a complex eigenvalue 
qe *« . We see that the time evolution of the coherent state \q > remains 
within the family of coherent states \q >. The property embodied in Eq. (1271) 
is the temporal stability of coherent states |^ > under the action of H . 

Thus, we conclude that the new coherent states \q > satisfy the Klauder's 
criteria set (i) - (iii) for generalized coherent states. 

Finally, let us introduce an alternative notation for |^ > in terms of the 
real q and imaginary p parts of \q >, that is, \q >= (q + ip) / y2H. Then from 
Eq. (TJj$ we have 

> = ]piq >= £ f^VV > • (28) 

The adjoint coherent states are defined as 



< C | =<p ,,| = £ < ^ m^pL f £±£T)>. (29) 
^ (2ft) n / 2 Vn! \ / 

Despite the fact that the adjoint state is labelled by q, the series expansion 
Eq. (1291) are formed in fact of powers of <;*. 

3.1 Quantum mechanical vector and operator repre- 
sentations based on coherent states \q > 

In spirit of Klauder's consideration [17], let's show that the resolution of unity 
criteria Eq. (fl9]) with W^(x) given by Eq. (125]) allows us to list fundamental 
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quantum mechanical statements pertaining to the associated representation 
of Hilbert space. Indeed, it is easy to see that the newly introduced coherent 
states |^ > provide: 

1. Inner Product of quantum mechanical vectors \tp > and \ip > defined 

as 

< y# >=\J d ^ < v?k > w^(kl 2 ) < >. ( 3 °) 

c 

where the vector representatives are wave functions < ip\q > and < q\ijj > 
given by 

< >= E ^^(^(-/^kl^)) 172 < ¥>|n >, (31) 

n=0 v 

< M >= E < n|V> > ^ T L (^ ) (-/^k| 2 ' 1 )) 1/2 - (32) 

n=0 

2. Vectors Transformation Law 

< q\A\*l> >=\J rf V < <^k' > ^(k'l 2 ) < ?V >, (33) 

c 

where < ?|»4|<r' > is the matrix element of quantum mechanical operator 

A. 

3. Operator Transformation Law 



< q\AiA 2 \q >=\ J d ^" < c|Ak" > ^(k'T) < ^"l^k" >, (34) 

c 

where A\ and A2 are two quantum mechanical operators. 

Further, the inverse map from the functional Hilbert space representation 
of coherent states \q > to the abstract one is provided by the following 
decomposition laws: 

4. Vector Decomposition Law 

IV>>=^/dVk> w^(kl 2 ) <<#> ■ (35) 

c 
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5. Operator Decomposition Law 



A= ~j d\ l( i\ 2 ■ \q x > W^tl 2 ) < qt\A\q 2 > W^ 2 \ 2 ) < ?2 |. (36) 

c 

Thus, we conclude that the resolution of unity Eq. ffl~9l) with W n (x), given 
by Eq. fl25l) . provides an appropriate inner product Eq. fl30l) and lets us intro- 
duce the Hilbert space. 

All of the above listed results lead to the well-know fundamental equations 
of quantum optics and coherent states theory [TTJ , [H] in the limit case /i = 1. 

Table 2 summarizes the definitions and equations attributed to the newly 
introduced coherent states \q > with those for the coherent states \z >. Table 
2 presents two sets of equations for a coherent state |... >, for an adjoint 
coherent state <...|, for the probability P(n) that n photons will be found in 
the coherent state |... >, for a positive weight function W(x) in the resolution 
of unity equations (IT91) and (1261) . the mean number < ...|a + a|... > of photons 
when the state is a coherent state |... >, and the quantum mechanical vector 
\ip > decomposition law. 





k> (o <n< l) 


\z > 


... > 


n=0 


oo 

n= ^ 


<... 


g<n|^r(4 B) (-^)) 1/2 




P(n) 




\*\ m e -\z\* 
nl 


W(x) 


1-M 1 

(1 — fi) • exp{ — (1 — n) *x} 


1 


< ... \a + a ... > 


(A*kr)/r(A* + i) 


kl 2 


|V > 


± f d\\q > W,W) < q\ip > 
c 


i / d 2 z\z >< z\ip > 
n c 



Table 2. Coherent states k > vs coherent states \z >. 



14 



4 Generalized Bell and Stirling Numbers 

4.1 Fractional Bell polynomials and Bell numbers 

Based on the fractional Poisson probability distribution Eq.([7j) we introduce 
a new generalization of the Bell polynomials 

= E""^^ a r(„(t + ») + !» ■ B -^°) = 1 - 

n=0 fc=0 VrV 'I'll 

(37) 

where the parameter /j is < /x < 1. We will call B^x, m) as the fractional 
Bell polynomials. 

The polynomials B^x, m) are related to the well-known Bell polynomials 
[TT] B(x, m) by 



x x r ' 



B^x,m)\, =1 = B(x,m) = e -*£ n m — . (38) 



n=0 



From Eq.( !37l) we come to a new formula for the numbers B^m), which 
we call the fractional Bell numbers 



^n m ^(k + n)\ (-l) k , . 

W = »A-)U = £ 7 £ l h r W + n) + i)) (39) 

n=0 fc=0 v ^ v 7 yy 

or 

oo m 

S M (m) = X)(-l) B ^ r 4 l) (-l), ( 4 0) 

n=0 

where Ej?\—x) is defined by Eq.([6]). The new formula Eq. (140p is in 
fact a fractional generalization of the so-called Dobinski relation known since 
1877 [15], [16]. Indeed, at n — 1 when the Mittag-Leffler function is just 
the exponential function, E\{x) = exp(— x), we have, (— l) n E± {— 1) = 
and the new equation (1401) becomes the Dobinski relation [15] for the Bell 
numbers B(m), 



B(m) = B,(m)\, =1 =e- 1 J2 1 T - ^ 



n! 

n=0 
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As an example, here are a few fractional Bell numbers 



B„(0) = 1, = TT7~~ 7T) £„(2 



1 _ , , 2 1 



r> + i)' ^ ' r(2// + i) r(/i + i) : 

n/n , 6 6 1 

5 m(3) = T7oT~T~T\ + F7TT-^T + 



r(3// + i) r(2// + i) r(/i + i)' 
24 36 14 1 



r(4/x + i) r(3/i + i) r(2/i + i) r(// + i)' 

Now we focus on the new general definitions given by Eqs. (l3"T|) and (HOI) to 
find the generating functions of the polynomials B^(x,m) and the numbers 
B^{m). Let us introduce the generating function of the new polynomials 



B^(x, m) as 



oo 

S 



F^x) = ^2— [ B^x,m). (42) 

m=0 

Therefore, to get the polynomial B^(x, m) we should differentiate F^(s, x) 
m times with respect to s, and then let s = 0. That is, 

gm 

B^(x,m) = —F^s,x)\ s=0 . (43) 

To find an explicit equation for F^(s,x), let's substitute Eq. fl37|) into 
Eq.( f42l) and evaluate the sum over m. As a result we have 

p , sn y^ (k + n)\ (-x) k 

W^-2^^[ e 2^ k\ T(fi(k + n) + l)Y 1 ] 

n=0 fc=0 V ^ V 1 11 

Then, introducing the new summation variable I — k + n, yields 

" ( ] h"' fc('-»)"v+i) 

2-<T^l + l)^n\{l-n)\ 1 J -Z.r^+1)- 
Finally, we obtain 
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F,(s,x) = E,(x(e s -l)), (45) 

where E^z) is the Mittag-Leffler function given by the power series 
Eq.©. 

It is easy to see that the generating function F^SjX), given by Eq. (j45|) . 
can be considered as the moment generating function of the fractional Poisson 
probability distribution (see Eq.(35) in Ref. [5]). 

In the case of fi — 1, Eq. (j4"5l) turns into the equation for the generating 
function of the Bell polynomials, 

F x (s, x) = exp{(a;(e s - 1))}. (46) 

If we put x = 1 in Eq. (l45p . then we immediately come to the generating 
function B^(s) of the fractional Bell numbers B^m) 

co m 

B ^ = E h B ^ m) = E ^ s ~ l) - (47) 

m=0 

The numbers B^m) can be obtained by differentiating B^{s) m times 
with respect to s, and then letting s = 0, 

B fl (m) = -^B^s,x)\ s=0 . (48) 
When ij,= 1, Eq.(jlZJ) reads 

00 m 

B ^ = E —M m ) = ex P{( e * - !)}» (49) 

m=0 

and we come to the well-known equation for the Bell numbers generating 
function. 

Let us show how the fractional Bell polynomials are related to the new 
coherent states |<j > introduced by Eq.f fl~4l) . For the boson creation a + and a 
annihilation operators of a photon field that satisfy the commutation relation 
[a, a + ] = aa + — a + a = 1, the diagonal matrix element of the n-th power of 
the number operator (a + a) n yields the fractional Bell polynomials of order 
n, 

<^\{a + a) n \q>= B^(\q\ 2 ,n). (50) 
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Then, the diagonal coherent state \q > matrix element < q\ exp(—iHt/h)\q > 
of the time evolution operator 

exp(—iHt/h) = exp {(— iuot / h)a + a} (51) 
can be written as 

oo .. 

< q\ exp {(-iut/h)a + a} \q >= £ -.(-i^)" < ^|(a+a) n |^ >= (52) 

ra=0 

n=0 

where /i is Planck's constant. 

By comparing with Eq. (j42p we conclude that 

< q\ exp {{-iuot/K)a + a} \q >= E^\q\ 2 (exp(-iujt/h) - 1)). (53) 

In other words, the diagonal coherent state \q > matrix element of the 
operator exp {(— iut / K)a + a] is the generating function of the fractional Bell 
polynomials. In the special case // = 1, Eq. (l53|) reads 

< z\ exp { (— iut/h)a + a} \z >= exp{(|^| 2 (exp(— iut/h) — 1)}, (54) 

that is, the diagonal coherent state \z > matrix element of the operator 
exp {(— iujt/K)a + a] is the generating function of the Bell polynomials. This 
statement immediately follows from Eqs. (11.5-2) and (11.2-10) of Ref . [1"9] for 
the diagonal matrix element of the operator exp {(— iuot / K)a + a} in the basis 
of the coherent states \z>. 

4.2 Fractional Stirling numbers of the second kind 

Now we are set up to introduce the fractional generalization of the Stirling 
numbers^] of the second kind S^(m, I) by means of equation 

4 Stirling numbers, introduced by J. Stirling [12j in 1730, have been studied in the 
past by many celebrated mathematicians. Among them are Euler, Lagrange, Laplace 
and Cauchy. Stirling numbers play an important role in combinatorics, number theory, 
probability and statistics. There are two common sets of Stirling numbers, they are so- 
called Stirling numbers of the first kind and Stirling numbers of the second kind (for 
details, see Chapter 8 in Ref. [T3] and [Tl]). 
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Bp(x, m) = y] S^m, l)x\ (55) 
1=0 

where B^x, m) is a fractional generalization of the Bell polynomials given 
by Eq. (!37|) and the parameter /i is < \i < 1. At // = 1, Eq. (l55l) defines 
the integers S(m,l) = S^m, i)U=ij which are called Stirling numbers of 
the second kind. At x = 1, when the fractional Bell polynomials B^x, m) 
become the fractional Bell numbers, B^{m) = B^x, m)\ x =i, Eq. (155l) gives 
us a new equation to express fractional Bell numbers in terms of fractional 
Stirling numbers of the second kind 

m 

B ll (m) = Y,S f +(m,l). (56) 

1=0 

To find S^m, I) we transform the right hand side of Eq. (!371) as follows 



,n ><■ r(A*(*+n) + l)) 

^ n! ^ 1 ; (/-n)!I>/ + l)' 1 ; 

here 9(1) is the Heaviside step function, 

Then, interchanging the order of summations in Eq.( l57j) yields 



OO 



L <<< ' ! = E T^TT — 77 E ^^777^7 = E Tv^T E(-^^ ' 
^T(ul + 1) t-^r n\(l-n)\ ^ T(ul + 1) 



Z=0 ~ vr "~ ' n=0 1=0 ~ vr " ' n=0 



(59) 

where the notation ( l ) = ,,/ ! has been introduced. 

\nJ n\(l—n)\ 

By comparing Eq. (1551) and Eq. (!59|) we conclude that the fractional Stirling 
numbers S^mJ) are given by 



^KO = f 7-7 TT ) E(-!) 1 - K> ( 6 °) 

" ' n=0 
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Sn(m, 0) = Smfl, S^(m,l) = 0, l = m + l, m + 2, .... 

As an example, Table 3 presents a few of fractional Stirling numbers of 
the second kind. 



m\l 


1 


2 


3 


4 


5 


6 


7 


1 


l 














2 


r(„+D 


r(2^+i) 












3 




6 

r(2/u+i) 


r(3/i+i) 










4 




14 




24 










r(2/Lt+i) 


r(3/i+i) 


r(4ju+i) 


5 


l 


30 


150 


24U 


120 








r(2^+i) 


r(3/H-i) 


r(4«+i) 


r(5At+i) 


6 




62 


540 


1560 


1800 


720 




ron-i) 


r(2M+D 


r(3M+i) 


r(4/x+i) 


r(5M+i) 




7 




126 


1806 


8400 




15120 


5040 


r(A»+i) 


r(2//+i) 


r(3/i+i) 


r(4ju+i) 


r(5^+i) 


r(6/i+i) 


r(7M+i) 



Table 3. Fractional Stirling numbers of the second kind S^mJ) (0 < 
^<1). 

Some special cases are 

S >' 1) = I>TT>' ^'"■ 2) = (2 ""- 1) f(2^Ty (61) 

S> ,3) = (3"-3.2" + 3 )f ^, (62) 
S„(n, 4) = (4» - 4 ■ 3" + 6 ■ 2" - 4) (63) 

s >-''- 1 »- 2 r('(l' -'% + !) ■ 5 * (,! ' n) = rO^W (64) 

It is easy to see that at \i — 1 the above equation turns into the well 
known representation for the standard Stirling numbers S(m,l) = Si(m,l) 
of the second kind [27] . 

n=0 ^ ' 
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Thus, one can conclude that there is a relationship between fractional 
Stirling numbers S^im, I) of the second kind and standard Stirling numbers 
S(m, I) of the second kind 

S»(m,l)= " + S(m,l). (65) 

or 

r(id + 1) 

/! 

Let's note that Eqs. fl65|) or fl66l) allow us to find new equations and iden- 
tities for the fractional Stirling numbers S^(m, I) based on the well-know 
equations and identities for the standard Stirling numbers S(m,l). 

To find a generating function of the fractional Stirling numbers S^m, I), 
let's expand the generating function F^(s,x) given by Eq.( H2i) . Upon substi- 
tuting B^x, m) from Eq. (1551) we have the following chain of transformations 



S( m , I) = »S M (m, I). (66) 



ml . 

m=0 \ 1=0 



OO / oo r 
S 



E-AE 9 (™- o^k i> 1 =EE s ^ ^ K> (67) 

m=0 \/=0 / 1=0 \m=l J 

where 6{m — I) is the Heaviside step function defined by Eq. (l58l) . 
On the other hand, from Eq. fl45p . we have for F M (s,x) 

Upon comparing this equation and Eq. (1^7|) , we conclude that 

E^K')^ = ^ T i ? ' = 0,1,2,.... (69) 

Now we are set up to define two generating functions Q^s, I) and J-^s, t) 
of the fractional Stirling numbers of the second kind, 
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m=l v ' 



(70) 



oo m m ,j oo s _ 1 V 

*) = EE s ^ o V = E tLtit = ^ (t(eS - (71) 

m=0 1=0 1=0 ^ ' 

As a special case /i = 1, equations flTOl and (ITTj) include the well-know 
generating function equations for the standard Stirling numbers of the second 
kind S(m,l) (for instance, see Eqs.(2.17) and (2.18) in Ref.[H]), 

Gi(s,l) = g„(s,l) | pF r%l)^ = {e ~ j , Z = 0,1,2,.... (72) 

*— ' m! /! 

and 



00 m m-d 

ftis, I) = ^(s, I) | M=1 = V V S(m, O^r = exp(t(e s - 1)), / = 0, 1, 2, .... 

z — ' z — 4 ml 

m=0 1=0 

(73) 

4.2.1 The Bernoulli numbers and fractional Stirling numbers of 
the second kind 

The Bernoulli numbers B n , n = 0, 1, 2, ...have the generating function [28] . 

E S A = ^T ( 74 ) 
' n! e { — 1 

n=0 

These numbers play an important role in the number theory. 
Let's show that the Bernoulli numbers B n can be presented in terms of 
the fractional Stirling numbers of the second kind S^(n,k) in the following 

way 

5n = E(-l) fc r(^ + l)^i. (75) 
fc=0 
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To prove Eg. (1751) we substitute B n from Eg. (1751) into the left hand side 
of Eg. (|74p . Therefore, we have 

f>^ = EEH)^+1)^ = (76) 
^— i n\ i f— i k + 1 n! 

^^(-i) fc r(^ + i)-^ (n ' fc) r 



n=0 n=0 fc=0 

oo oo 



k + 1 n\ 

k=0 n=k 



To transform the right hand side of Eg. (l76l) we use Eg.( j70l) and obtain 



OO 



y bJ-- = y y(-i) fc r(^+i)^^l^ = y^if^zli 

^ n\ 1 VP ' k + 1 n\ ^ K ' k + 1 e'-l 

n=0 k=0 n=k k=0 

Thus, we have proved Eg. (1751) . 

In the case \i = 1, when the gamma function is T( / uA;+l)| At=1 = k\, Eg. (1751) 
reads 

k=0 

and we recover the representation of the Bernoulli numbers B n in terms of 
the Stirling numbers of the second kind S(n, k) (for instance, see the eguation 
for B n on page 2547 of Ref. [1] 



5 Statistics of the fractional Poisson proba- 
bility distribution 



Now we use the fractional Stirling numbers introduced by Eg. (l60|) to get 
the moments and the central moments of the fractional Poisson process with 
the probability distribution function P M (n, t) given by Eg.([7j). Indeed, by 
definition of the m-th order moment of the fractional Poisson probability 
distribution we have 



n 



oo 
n=0 



n 



n=0 



nl 



E 

k=0 



k\ T(fi(k + n) + l)Y 



< n < 1. 
(78) 
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It is easy to see that n™ is in fact the fractional Bell polynomial B^ut^, m) 
defined by Eq.(j57j). 

From other side, with the help of Eqs. fl37|) and fl55|) we find 



11" 



J2S,(m,l)(utn l - (79) 



m 

,m 

Z=0 

Hence, the fractional Stirling numbers S* M (m, Z) of the second kind natu- 
rally appear in the power series of the m-th order moment of the fractional 
Poisson probability distribution. 

Using analytical expressions given by Eqs. fjoDl) and (IT9"j) . let's list a few 
moments of the fractional Poisson probability distribution 

% = X> P >'*) = rfc + i) ' (80) 



n=0 



n 2 = \^ n 2p , / t = v — ; + 81 

m , ; r(2/x + i) r(/x + i)' v ' 



n=0 

00 



,3 



n- = > n 3 P u n, t) = — ^ ^- + — ^ ^- + — r, 82 

" r(3// + i) r(2// + i) r(// + i)' 1 ; 



S^^t)^^^^. (83) 
" ^ r(4/x + i) r(3 M + i) r(2 M + i) r^ + i) 1 ; 

In terms of the power series over the first order moment n^, the above 
equations (IHTI) - (!83|) read 



"* = r(2„ + i) + (84) 

eyq. + iy 6(ry + iy 
r(3/< + i) " + r(2 M + i) "» + "<" lS5) 

= 24 ( r ( „ + 1)) V sqrfr + lyf i4 ( r ( „ + 1£ 

r(4/x + i) r(3/i + i) r(2/x + i) ^ " v 7 
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The mean n M Eg. (I80|) and the second order moment were first obtained 
by Laskirjf] (see, Eqs.(26) and (27) in Ref.0). 

In the case when jj, — 1, equations - fl86|) become the well-know 
equations for moments of the standard Poisson probability distribution with 
the parameter n = ri\ = vt (for instance, see Eqs.(22) - (24) in Ref. [29J ) . 

5.1 Variance, skewness and kurtosis of the fractional 
Poisson probability distribution 

To find analytical expressions for variance, skewness and kurtosis of the frac- 
tional Poisson probability distribution, let's introduce the central m-th order 
moment M M (m) 

oo 

M M (m) = (rv - n„) m = J^(n - n^P^n, t) = 

n=0 

oo m / 

ED- 1 )" u r^ rp ' ,(M) = (87) 

n=0 r=0 ^ ' 

m / \ r 

E(- 1 ) m_p 7 { ^ r ~ r S ^ 

where S^(r, I) is given by Eq.( l60l) . 

Hence, in terms of power series over the first order moment Eq. (1801) . we 
have 

M„(l) = 0, (88) 
M(3) 2 p( r ^ +1 » 3 3 ' r ^ +1 >' 2 |lW (90) 

r(3/1 + 1) r(2/i + i) +1 j " 1 



5 The second order moment defined by Eq. (f8Tj) can be presented as Eq.(27) of Ref. [5] if 
we take into account the well-know equations for the gamma function T(/i) 

r( M + 1) = fjT(fi), r(2 M ) = —=-t(jm) • r(p + -). 

\/7T 2 
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2(r(// + i)) 



2 



r(2// + i 



3 1 '. 7/ -1 lnf t + n 



> i '"fii 



M (4) = 3 ( 3(r(/i + 1))4 - 8(F(/i + 1))3 + 4(IV + 1))2 - 1 ) n*+ (91) 

Mi J dl r(4/i + i) r(3// + i) + r(2/z + i) 1 » + " ' 



6 /6(T(m + 1))' _ 4(r(, + !))» \ /7(T(m + 1))' _ N 

V r(3// + i) r(2/x + i) / \ r(2/i + i) / " " 

Further, in terms of the above defined central moments M M (m), the vari- 
ance a 2 , skewness s^, and kurtosis of the fractional Poisson probability 
distribution are 

o* = M,(2), (92) 
M M (3) 

— o/n ; ~j (93) 
M^ /2 (2) 

fc = —£-1 - 3. 94 

M2(2) v 7 

In the case when /i = 1, new equations fl89|) - fl9TT) turn into equations for 
the central moments of the standard Poisson probability distribution (see, 
Eqs.(25) - (27) in Ref. |29j). 

Equations f[9~2"j) - ([BP , at /i = 1, turn into the equations for variance, 
skewness, and kurtosis of the standard Poisson probability distribution with 
the parameter n (for instance, see Eqs.(29) - (31) in Ref. |29j). 

Table 4 summarizes equations for the Bell polynomials B(x,m), the Bell 
numbers B(m), generating function of the Bell numbers B(s), the Stirling 
numbers of the second kind S(m, I), generating function of the Stirling num- 

oo 

bers of the second kind ^2 S^mJ)^, and the m-th order moment, at- 

m=l 

tributed to the fractional Poisson distribution with those for the standard 
Poisson distribution. 
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fractional Poisson (0 < /i < 1) 


Poisson (// = 1) 


B(x, m) 


V ^mz™ V (k+n)\ (-x) k 

2^ 11 „! Z-/ fc! r0i(fc+n)+l)) 
n=0 fc=0 


-x y 
n=0 


B{m) 


n=0 


n=0 


Bis) 


^(e s - 1) 


exp{(e s — 1)} 


S(m, I) 


rwi)E(-iy-"Ci)^ 

n=0 


£E(-i)<- n O m 

n=0 


oo 
m=l 


(e s -l)< 

r( M /+i) 


(e s -l)< 

n 




«=0 


1=0 



Table 4. Polynomials, numbers, moments and generating functions at- 
tributed to the fractional Poisson process vs the standard Poisson process. 



6 Conclusions 

The quantum physics and the number theory applications of the fractional 
Poisson probability distribution have been developed. 

The key new results related to quantum physics applications are given by 

Eqs.flU, (USD, PU), (ESD> dSU) - H35D- 

The key new results related to the number theory applications are given 

by Eq S .(j37D, m, (USD, (33), dSSD, (HQ}, ([HD, (USD, CSD, dHOD - dH, 
(EZD - JSU). 

Tables 2,3,4 summarize our findings for the fractional Poisson probability 
distribution in comparison to the well-known results related to the standard 
Poisson probability distribution. 

These findings facilitate the further exploration of other fields to apply 
the fractional Poisson probability distribution and will provide deep insight 
into long-memory impacts on counting processes. 
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7 Appendix 



To obtain Eg. (1241) we use the Laplace transform of the Mittag-Leffler function 

oo 

dre-TE^-ZT*) = -i-, 
i- ~\~ z 



for instance, see equation (26) on the page 210 in Ref . [25] . 
Changing the variable r — > st and the parameter zs^ — > ( yields 



oo 

/■ 



S 



dte-°%(-(t») = —. (95) 



o 



By differentiating Eq.( l95i) n times with respect to ( we obtain, 

oo 

dte~^E^\-^)= ^ + c)n+l . (96) 

o 

It is easy to see that at ( = \i Eq.(|96l) goes into Eq.( !24l . 



References 

[1] N. Laskin, Fractional quantum mechanics and Levy path integrals, Phys. 
Lett. A268 (2000) 298-305. 

[2] N. Laskin, Fractional quantum mechanics, Phys. Rev. E62 (2000) 3135- 
3145 (available on-line at http://arxiv.org/abs/0811.1769). 

[3] R.P. Feynman, The Space-Time Formulation of Nonrelativistic Quan- 
tum Mechanics, Rev. Mod. Phys. 20 (1948) 367-387. 

[4] R. P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals, 
McGraw-Hill, New York, 1965. 

[5] N. Laskin, Fractional Poisson process, Communications in Nonlinear 
Science and Numerical Simulation, 8 (2003) 201-213. 



28 



M. Kuno, D.P. Fromm, H. F. Hamann, A. Gallagher, and D.J. Nesbitt, 
Nonexponential "blinking" kinetics of single CclSe quantum dots: A 
universal power law behavior, J. Chem. Phys. 112 (2000) 3117-3120. 

L. Sabatelli et al, Waiting time distributions in financial markets , The 
European Physical Journal B - Condensed Matter and Complex Sys- 
tems, 27 (2002) 273-275. 

W. Willinger, and V. Paxson, Where Mathematics Meets the Internet, 
Notices of the American Mathematical Society, 45 (1998) 961-970. 

CM. Zaslavsky and A. I. Saichev, Fractional kinetic equations: solutions 
and applications, Chaos 7 (1997) 753-764. 

J.R. Klauder, The Current State of Coherent States, 



arXiv:quant-ph/0110108fr rl 17 Oct 2001. 

E. T. Bell, Exponential Polynomials, The Annals of Mathematics, Sec- 
ond Series, 35 (1934), 258-277. 

J. Stirling, Methodus Differentialis: Sive Tractatus de Summatione et 
Interpolatione Serierum Infinitarum, Gul. Bowyer, London, 1730. (En- 
glish translation by J. Holliday, The Differential Method: A Treatise of 
the Summation and Interpolation of Infinite Series, 1749). 

C. A. Charalambides, Enumerative combinatorics, Ch.8, Chap- 
man&Hall/CRC, 2002. 

C. A. Charalambides and J. Singh: A review of the Stirling numbers, 
their generalizations and statistical applications, Communications in 
Statistics - Theory and Methods, 17 (1988) 2533-2595. 

G. Dobinski, Summierung der Reihe ^n m /n!, fiir m = 1,2,3,4,5,..., 
Grunert Archiv Arch, fiir M. und Physik 61 (1877) 333-336. 

Gian-Carlo Rota, The Number of Partitions of a Set, American Mathe- 
matical Monthly, 71 (1964) 498-504. 

J.R. Klauder, E.C.G. Sudarshan, Fundamentals of Quantum Optics, 
Benjamin, New York, 1968. 



29 



R. J. Glauber, Quantum Theory of Optical Coherence. Selected Papers 
and Lectures, Wiley- VCH, Weinheim 2007. 

L. Mandel and E. Wolf, Optical Coherence and Quantum Optics, Cam- 
bridge University Press, 1995. 

K.B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, 
New York, 1974. 

S.G. Samko, A. A. Kilbas, and O.I. Marichev, Fractional Integrals and 
Derivatives and Their Applications, Gordon and Breach Science Pub- 
lishers, Langhorne, PA, 1993. 

K.S. Miller, An Introduction to the Fractional Calculus and Fractional 
Differential Equations J. Wiley & Sons, New York, 1993. 

C.W.Gardiner, Handbook of Stochastic Methods, Second Edition, 
Springer Series in Synergetics, Editor H Haken, 13, Springer- Verlag, 
1985. 

G. Mittag-Leffler, Sur la representation analytique d'une branche uni- 
forme d'une fonction monogene, Acta Mathematics, 29 (1905) 101-182. 

A. Erdelyi, Ed., Higher Transcendental Functions, Vol.3, pp. 206-227 
(Chapter 18 Miscellaneous Functions) McGraw-Hill, New York, 1955. 

E. Schrodinger, Der stetige Ubergang von der Mikro-zur 
Makromechanik, Die Naturwissenschaften 14 (1926) 664-666. 

M. Abramowitz and LA. Stegun, (Eds.). Stirling Numbers of the Second 
Kind. §24.1.4 in Handbook of Mathematical Functions with Formulas, 
Graphs, and Mathematical Tables, 9th printing, pp. 824-825, New York, 
Dover, 1972. 

M. Abramowitz and LA. Stegun, (Eds.). Bernoulli and Euler Polynomi- 
als and the Euler-Maclarin Formula §23.1 in Handbook of Mathematical 
Functions with Formulas, Graphs, and Mathematical Tables, 9th print- 
ing, pp. 804, New York, Dover, 1972. 

EW. Weisstein, Poisson Distribution, From MathWorld - A Wolfram 
Web Resource, 



http: //mathworld. wolfram.com/PoissonDistribution. html 



30 



